Thermally induced instability of a doubly quantized vortex in a Bose-Einstein 

condensate 
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We study the instability of a doubly quantized vortex topologically imprinted on 23 Na condensate, 
as reported in recent experiment [Phys. Rev. Lett. 93, 160406 (2004)]. We have performed nu- 
merical simulations using three-dimensional Gross-Pitaevskii equation with classical thermal noise. 
Splitting of a doubly quantized vortex turns out to be a process that is very sensitive to the presence 
of thermal atoms. We observe that even very small thermal fluctuations, corresponding to 10 to 15% 
of thermal atoms, cause the decay of doubly quantized vortex into two singly quantized vortices in 
tens of milliseconds. As in the experiment, the lifetime of doubly quantized vortex is a monotonic 
function of the interaction strength. 
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Experimental studies of vortices in Bose-Einstein con- 
densates have revealed their peculiar properties related 
to the quantized circulation originally predicted by 
Onsager and Feynman in the context of rotating super- 
fluid 4 He and further explored in Refs. Q and Q. 
The quantization of the circulation is a manifestation of 
the existence of a macroscopic wave function. Direct 2ir- 
change of the phase of the condensate wave function when 
going around the vortex core was experimentally demon- 
strated by using the interferometric technique @. The 
quantized vortex can not just disappear, it can leave the 
condensate or annihilate with a vortex having the oppo- 
site circulation. There exists another route for vortices 
that have multiple topological charge, they can also split 
into several vortices having smaller charges decreasing in 
such a way the energy of the system. 

In recent experiments 0, Q doubly quantized vortices 
have been imprinted in Bose-Einstein condensate of 23 Na 
atoms. A novel approach allowing to create vortices 
with multiple topological charge was implemented @. 
In this new method, as opposed to dynamical phase- 
imprinting techniques like rotating the atomic cloud in 
an anisotropic trap or stirring the condensate with the 
help of a laser beam, vortices are generated by adiabat- 
ically reversing the magnetic bias field along the trap 
axis. This topological phase imprinting technique leads 
to vortices displaying winding numbers 2 or 4 depending 
on the hyperfine state the sodium condensate was pre- 
pared in (|1, —1 > and |2, +2 >, respectively) [g. In Ref. 
7] the authors investigate the evolution of doubly quan- 
tized vortices by using a tomographic imaging method 
@ . They observe the decay of doubly quantized vortex 
into singly quantized vortices and suggest the possible 
explanation of this splitting as being a result of dynam- 
ical instability, however, not specifying the character of 
the perturbation seeding the instability. 

Existing theoretical explanation of decay of doubly 
quantized vortices involves the analysis of stability of 
the vortex at zero temperature in terms of cigenmode 
spectrum of the Bogoliubov equations in two- and three- 
dimensional geometry as well as the numerical solution of 
the Gross-Pitaevskii equation in slightly anisotropic trap 



|lfj| . Two-dimensional Bogoliubov eigenvalue spectrum 
shows complex f requ encies for certain values of the inter- 
action strength |ToL ITlT|. In fact, quasiperiodic behavior 
is found - the stability windows are followed by the in- 
stability regions. Neither the anisotropy of the trap nor 
the phenomenologically introduced dissipation is able to 
force the decay of doubly quantized vortex when there are 
no complex eigenvalues (i.e., the stability window condi- 
tions are fulfilled) |l0j . However, in the experiment the 
monotonic increase of the lifetime of the vortex with the 
strength of the nonlinearity is observed Q ■ It turns out 
that such behavior is attributed to the three-dimensional 
geometry of the system and its origin could be due to the 
presence of trap anisotropy (as suggested in Ref. ^(|) or 
the thermal noise (as claimed by this Letter). 

In this Letter we show that only thermal fluctuations, 
not other disturbances (as, for example, due to the con- 
finement anisotropy), lead to the decay times compara- 
ble with that reported in the experiment Q. To this 
end, we have performed numerical simulations using the 
three-dimensional Gross-Pitaevskii equation in the ver- 
sion described in Ref. [l^, i.e., with a classical thermal 
noise. We find that even extremely small thermal fluctu- 
ations dramatically accelerate the decay of doubly quan- 
tized vortex into two singly quantized vortices. Increas- 
ing the number of thermal (uncondensed) atoms already 
to 10 to 15% reduces the lifetime of the vortex below 
100 ms. Therefore, we argue that although the authors 
of Ref. report that the experiment is performed under 
the condition of no discernible thermal atoms presence, 
the decay of the vortex is triggered by thermal rather 
than quantum fluctuations. 

We describe the system of degenerate bosonic atoms 
in terms of the classical field that is the complex func- 
tion representing all atoms not only those within the 
condensed faction. At zero temperature all atoms are 
condensed and the classical field becomes the conden- 
sate wave function which satisfies the Gross-Pitaevskii 
equation. However, as it was argued in Ref. |12|. the 
same equation is fulfilled by the classical field at finite 
temperatures, although in this case it must be inter- 
preted in a different way. The high energy solution of the 
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Gross-Pitaevskii equation describes both the condensed 
and thermal atoms. To get correct physical interpreta- 
tion of the classical field one has to perform averaging 
over time or space of corresponding single-particle den- 
sity matrix. The condensate wave function is just the 
dominantly populated eigenmode of the averaged single- 
particle density matrix. Other modes represent thermal 
atoms. 

We start our simulations with the wave function of the 
cigar-shaped condensate with a vortex imprinted along 
the symmetry axis. The topological charge of the vor- 
tex equals 2. The trap parameters are the same as in 
the experiment of Ref. 7], i.e., the radial (axial) trap 
frequency equals 220 Hz (3 Hz). The number of atoms is 
changed in such a way that the values of the control pa- 
rameter defined as an z , where a is the scattering length 
and n z = J | ip(x, y, z — 0) \ 2 dx dy is the linear atom den- 
sity along the symmetry axis taken at the center of the 
trap, varies up to about 14 (see Ref. Q). To find the 
condensate wave function with the vortex we solve the 
Gross-Pitaevskii equation in imaginary time. The nu- 
merical code im plem ents an operator splitting technique 
as described in |l3| and works on two-dimensional grid 
since the phase of the vortex state, e l2< ^, can be separated 
off analytically. 

Next, we randomly disturb the condensate wave func- 
tion, i.e., introduce the thermal noise into the classical 
field. Since the classical field is a complex function we 
perturbed independently the amplitude and the phase. 
Both approaches give very similar results. The evolution 
of the classical field ip is governed by the time-dependent 
Gross-Pitaevskii equation (see Ref. [l2T |) 
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V 2 + V trap +gN\iP\ 2 )iP, (1) 



where g = 4Trh 2 a/m. The above equation is solved on a 
three-dimensional grid using the Fast Fourier Transform 
split-operator method. We monitor the diagonal part of 
the single-particle density matrix that is just the density 
of the classical field. 

In Fig. we plot the time needed to split the dou- 
bly quantized vortex as a function of the energy pumped 
into the system while introducing the thermal noise. For 
each curve (i.e., the value of the interaction strength) 
the considered energies range from the zero temperature 
energy to an amount approximately twice larger. We 
have checked that this corresponds at most to 20% con- 
densate depletion. The basic observation is that the de- 
cay time is extremely sensitive to the presence of the 
noise. The rise of the energy by less than 1% (in com- 
parison with the zero temperature energy) already results 
in a violant decrease of the decay time. For example, it 
drops from 200 ms to 87 ms for an z — 5.9 (the case of 
zero temperature energy per particle equal to 250 ftw z ). 
Rigorously speaking, at zero temperature and without 
any disturbances that might initialize the decay (like the 
trap anisotropy or any deficiencies related to the process 
of phase imprinting) the lifetime of the vortex should 
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FIG. 1: Decay time of doubly quantized vortex as a function 
of the energy per particle. The successive curves correspond 
to different values of the control parameter defined as an z , 
where a is the scattering length and n z is an axial density at 
the center of the trap. From top to bottom an z equals 12.2, 
10.0, 7.2, 5.9, 3.9, and 2.3, respectively. 
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FIG. 2: Density cuts at zero temperature and an z = 10.0 
for different times showing that even at zero temperature 
the doubly quantized vortex eventually decays and two singly 
quantized vortices become disentangled along the whole con- 
densate. 



be infinite. Indeed, curves in Fig. ^ show well visible 
asymptotes when approaching zero temperature energies 
(although the top-lying end of each curve is a zero tem- 
perature decay time that is finite due to the presence 
of numerical noise). Further increase of the amount of 
the thermal noise causes additional, although slower, de- 
crease of the decay time. Fig. ^ shows that the decay 
time becomes comparable with the experimentally mea- 
sured values (tens of milliseconds) only when the thermal 
noise at the appropriate level is included in the dynamics. 

Even at zero temperature the vortex decays (see Fig. 
121 . The decay process sets in independently of the density 
as opposed to what is predicted based on two-dimensional 
spectral analysis of the Bogoliubov eigenmodes. Accord- 
ing to this analysis, the stability windows with respect 
to the interaction strength exist, the first one appear- 
ing approximately for 3 < an z < 12. They are defined by 
the requirement that all eigenfrequencies of all excitation 
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modes are real. Within the stability window no decay of 
the doubly quantized vortex is expected. Outside the 
window, i.e., when there is a complex eigenfrequency, 
the corresponding excitation mode grows exponentially 
in time resulting in the splitting of the vortex. Adding 
the third dimension changes the results dramatically as 
it is shown in Figs. Q an d Splitting of the vortex in 
three-dimensional condensate need not be a simple and 
uniform process. In fact, it starts in the region where the 
local value of an z is in the two-dimensional instability 
window (see Ref. [HJ)- So, the vortex splitting begins 
at the ends of the condensate (assuming an z < 12) and 
propagates towards the center of the trap. In Fig. 
we plot the lifetime of the doubly quantized vortex un- 
derstood just as the time needed to split the vortex and 
disentangle it along the whole condensate. 
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FIG. 3: Tomographic images of the condensate that illustrate 
the dependence of the decay time on the energy of the system. 
The interaction strength equals an z = 3.9 and the energy 
per particle is 234, 235, and 331 hui z from top to bottom, 
respectively. The lowest panel shows the separation of two 
singly quantized vortices in the case when E/N = 235fio; z 
(squares) and E/N = 331 hu z (triangles). 

In Fig. |3]we compare the decay time of a doubly quan- 
tized vortex for different levels of the introduced thermal 
noise. The pictures are obtained according to the tomo- 
graphic imaging technique as described in Ref. |7|. In 
this method, atoms within a 30 /im thick central slice 
of the condensate are transfered to a different hyperfine 
state and then imaged by using another resonant laser 
pulse. Following this prescription, each frame in Fig. |3| 



shows the contour plots of the density integrated axially 
within the slice of 30 /im thickness located at the center 
of the trap. The first two rows demonstrate how impor- 
tant is the thermal noise. Introducing the thermal noise 
on a very low level (the energy of the system is increased 
by less than 1%) already decreases the time needed to 
split the vortex by half. Rising the level of the thermal 
noise causes further lowering of the decay time what is 
necessary to get an agreement with the experiment. In 
the bottom panel we plot the separation between two 
singly quantized vortices as a function of time. The case 
when the energy of the system is increased by approx- 
imately 50% in comparison with the zero temperature 
energy is denoted by triangles. Although the population 
of the condensate is now about 94%, the distance between 
the two cores (~ 3 /im) almost does not change over the 
period of 90 ms. 
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FIG. 4: Decay time of doubly quantized vortex as a function 
of the density. Three sets of points correspond to different 
values of the energy brought in the system while introducing 
the thermal noise. From top to bottom the relative increase of 
the energy equals 40% (squares), 85% (triangles), and 135% 
(stars). The condensate depletion amounts to 6%, 13%, and 
20%, respectively. 

Finally, in Fig. ^ we present the data regarding the 
time of the decay of doubly quantized vortices which are 
intended to reproduce the main result of the experimental 
work • The authors of Ref. say only that the exper- 
iment was performed at the lowest possible temperature. 
Since the temperature of the system is determined based 
on the expansion of the thermal cloud, this statement 
puts on, in fact, the constraint on the number of thermal 
atoms. To our knowledge the presence of less than 15% 
of thermal atoms can not be detected by using currently 
available experimental techniques. Therefore, in Fig. 0] 
we plot three sets of data, each of them corresponding to 
different level of thermal noise. More precisely, the rela- 
tive increase of the energy of the system is constant for 
each data set. It is clear from Fig. 0]that the decay time 
gets larger when the interaction strength increases what 
remains in agreement with the experiment. For larger 
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density, however, the decay time is getting saturated. 
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FIG. 5: Density cuts along the radial plane z — at zero (the 
first row) and finite (the bottom row) temperatures for the 
interaction strength an z = 3.9. 

To answer the question whether other sources of in- 
stability are able to reproduce the experimental results 
of Ref. we investigated the influence of the trap 
anisotropy (as high as 2.3% according to 01) caused 
by the gravitational sag as well as imprinting deficiencies 
related to the rapid change of the magnetic field produc- 
ing the vortex and resulting in radial squeeze and vertical 
kick of the condensate 7. Il4| . We checked that the trap 
anisotropy makes the lifetime of a doubly quantized vor- 
tex finite but still bigger than the experimental values. 
For example, for an z — 10.0 the decay time equals 170 ms 
whereas for an z = 3.9 it is 76 ms. Moreover, since the 
gravitational field do not fluctuate the trap anisotropy 
can not explain the huge scatter of data in Fig. 3 of 
Ref. 0- Shifting the whole condensate off the center of 
the trap also does not help. It makes the total energy 
higher by increasing the potential energy not the kinetic 
one. Therefore the number of uncondensed atoms re- 
mains constant and the lifetime of the vortex equals that 
for unshifted condensate. Finally, we considered the dy- 
namical phase imprinting instead of the topological one. 
In this case the energy pumped into the system goes to 
the kinetic energy and results in a production of uncon- 



densed atoms. The decay time is comparable with the 
experimental values revealing in this way the basic role 
of the uncondensed atoms in the process of decaying of 
doubly quantized vortices. 

Another interesting phenomenon is discovered while 
tracing the dynamics of the splitting and subsequent evo- 
lution of two singly quantized vortices for longer times. 
It turns out that additional vortices enter the conden- 
sate. In Fig. [5] the first row corresponds to the case 
of zero temperature. Here, we observe that two singly 
quantized vortices (as numerically checked by using the 
interference technique as described in |l5j ) enter the con- 
densate at the same time and settle into a lattice (in a 
rotating frame of reference) realizing the first scenario re- 
ported in Ref. 0] . At finite temperatures (the last row) 
the dynamics is reacher. First, two additional vortices 
simultaneously enter the condensate, later on, however, 
one vortex is lost and the system ends with the lattice 
consisting of three vortices. 

In conclusion, we have addressed the issue concern- 
ing the lifetime of a doubly quantized vortex raised by 
a recent experiment 7j. We show that the decay of a 
doubly quantized vortex is driven rather by the thermal 
fluctuations than other kinds of perturbation (quantum 
fluctuations since we need 10 to 15% of thermal atoms 
to achieve the agreement with experiment |jj or the trap 
anisotropy). The thermal noise is the only reason, we 
find, able to decrease the time needed to split the vortex 
and make it comparable to the experimental values. As 
a result, the decay time is a monotonic function of the 
interaction strength. Our stressing of the role of thermal 
noise docs not contradict the statement of the authors of 
Ref. saying that the experiment was performed at the 
lowest possible temperature. 
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